/vj/vsA' ca - >m. *51^ 


NASA Contractor Report 189575 
ICASE Report No. 91-85 


ICASE 


CONVERGENCE OF CONTROLLERS DESIGNED 
USING STATE SPACE METHODS 


K. A. Morris 


Contract No. NAS 1-1 8605 
December 1991 

Institute for Computer Applications in Science and Engineering 
NASA Langley Research Center 
Hampton, Virginia 23665-5225 

Operated by the Universities Space Research Association 


National Aeronautics and 

Space Administration 

Langley Research Center 
Hampton, Virginia 23665-5225 




NAS A-CR- 189575 
19920011420 




Convergence of Controllers Designed Using State Space Methods 


K. A. Morris 1 

Department of Applied Mathematics 
University of Waterloo 
Waterloo, Ontario N2L 3G1 


ABSTRACT 

In this paper convergence of finite-dimensional controllers for infinite-dimensional sys- 
tems designed using approximations is examined. Stable coprime factorization theory is 
used to show that under the standard assumptions of uniform stabilizability/detectability, 
the controllers stabilize the original system for large enough model order. The controllers 
converge uniformly to an infinite-dimensional controller, as does the closed loop response. 
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1 Introduction 

We consider semigroup control systems on a Hilbert space X: 

x(t ) = Ax(t) + Bu(t ), (1) 

y(t) = Cx(t ), 

®(0) = x 0 , x 0 G D(A ) C X. 

The operator A generates a strongly continuous semigroup of operators 5(f) on X so A is 
closed with domain T>(A) dense in X. We assume bounded control and observation, and 
that the input and output spaces are finite-dimensional: B € B(R m ,X), C € B(X, R p ). This 
control system will often be abbreviated (A, B, C). For control functions in £ 2 ( 0 , 00 ; R m ) we 
interpret the solution of (1) in the mild sense: 

y(t) = CS(t)xo + C f S(t - T)Bu(r)dr. 

JO 

Further details can be found in, for instance [14]. 

The equations (1) are a model for a number of control problems, including those where 
the system dynamics are described by partial differential equations and hereditary differential 
systems. A closed form solution can be computed only in the simplest of situations. In 
general, it is necessary to use an numerical approximation to the semigroup control system 
(1) in order to simulate the response of the system . This approximation will typically be a 
system of n ordinary differential equations, which we write as 

x(t) = A n x(t) + B n u(t), (2) 

y(t) = C n x(t), 
x(0) = x 0n - 

This approximation is also used to compute controllers for the original system (1). Further 
details on the approximation scheme will be presented in a subsequent section. The above 
finite-dimensional system will often be abbreviated as ( A n , B n , C n ). 

In this paper, we are concerned with finite- dimensional controllers of the form 

z(t) - (A n - B n K n )z(t) + F n (y(t ) - C n z(t)) (3) 

Wn(t) = K n z(t ). 

There are a number of convergence questions associated with this approach: 
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• Does the controller (3) stabilize (1) for large enough n ? 

• Do the controllers (3) converge to an infinite- dimensional controller? 

• What relation does the closed loop performance of the controller (3) with (A n , B n , C n ) 
have to its implementation with the original system (A, B, C)1 

Convergence in the graph topology [18, 19] of the approximating systems to the origi- 
nal system is necessary for the validity of an approximation scheme as a basis for controller 
design. Convergence in the graph topology is equivalent to convergence of a sequence of 
coprime factors ( N n ,D n ) for the approximations to a coprime factorization for the original 
system. Without such convergence it cannot be concluded that a controller which stabi- 
lizes ( A,B,C ) also stabilizes (A„, B n ,C n ) for sufficiently large n, or that the closed loop 
responses converge. In fact, in the absence of convergence in the graph topology, the set of 
controllers which stabilize both the infinite-dimensional system ( A , B , C ) and the approxi- 
mations (A n , B n , C n ) may be empty, even for large n. 

In this paper this approach is used to obtain complete answers to the above questions 
for a wide class of systems. A complete answer to the problem of convergence of LQG type 
controllers is given. Although the graph topology is the fundamental topology underlying 
these convergence questions, it is not used explicitly in this paper. All results are derived 
directly using only theory on Bezout, or stable coprime factors. First, some background 
material on stable coprime factorizations is given. 

2 Stable Coprime Factorizations 

Suppose a system of the form (1) with Xo = 0 maps inputs in £ 2 ( 0 , 00, R m ) to outputs in 
£ 2 ( 0 , 00, R p ), and that furthermore, there is a maximum ratio, the £ 2 -gain between the norm 
of the output and the norm of the input: 

II ^ II £2(0, 00 , HP) — 7ll u ll,C 2 (0,oo,R m )- 

Then the system is said to be £ 2 -stable, and by the Paley- Weiner Theorem, the Laplace 
transform of such a system is a matrix with entries in . Here Tioo indicates the Hardy 
space of functions G(s) which are analytic in the right-half plane Re(s ) > 0 and for which 

sup sup \G(x +ju»)j < 00. 
w ®>o 

The norm of a function in Hoo is 

II G lloo = supsup|G(s+ jw)|. 

to x>0 
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We denote matrices with entries in 'H co by M( 7-foo ). The norm of a function in M( 'H< x > ) 
is the induced matrix norm 

II GIL = sup sup a [G{x + jut)] 

x>0 u 

where a denotes the largest singular value. 

The Laplace transform of a linear time-invariant system such as (1) is called its 
transfer function and is given by CR(s-, A)B where R(s\ A) indicates the resolvant of A. 
Suppose G is the transfer function of a given system, for which we wish to design a controller 
with transfer function H, of compatible dimensions, arranged in the feedback configuration 
shown in Figure 1. 
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Figure 1. Feedback System 


The 2x2 transfer matrix A((7, H) which maps the pair (ri, r 2 ) into the pair (ei, e 2 ) is given 

by 

' (I + GH)~ l —G(I + HG)~^ ‘ 

A (G,ff) = 

H(I + GH )-» (I + HCy 1 

The feedback system, or alternatively the pair ( G,H ), is said to be externally stable if each 
of the four elements in the above matrix belongs to the set S of stable transfer functions. We 
could define stability in terms of the transfer matrix from (r i, r 2 ) to (2/1 , 2 / 2 ); both notions of 
stability are equivalent [19]. Definition of S depends upon the application. Thus the closed 
loop system is ZVstable if and only if all four elements belong to M( 'H rx , ). The set of all 
plants which stabilize G is written S(G ): 

S(G) = {H | A (G,H) e M{ H 0 0 )}. 

Note that the present definition of stability is symmetric in G and H. Thus G stabilizes H 
if and only if H stabilizes G. 

For the common situation where the system is already stable and an aim of controller 
design is to improve the settling time of the system, we specify a real number cr > 0 which 
is the minimum acceptable stability margin. Then a system is said to be input/output 
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a -stable if its shifted transfer function G{s — a) is in M( Hoo ). Equivalently, we replace 
H 0 o by the algebra of functions which are analytic in the right half plane Re(s) > —a, 

and for which 

sup sup | G(x + ju>) | < oo 

U> 07>“‘<7 

with corresponding norm. 

Much of modern control theory is concerned with coprime factorizations of systems. 
The transfer function of a possibly unstable system G is written as the ratio of two coprime 
stable systems. For the case of £ 2 -stability, the transfer function of a system is written as 
G = ND -1 where N, D £ M( Hoo ) and there exists X, Y £ M( TCoo ) with 

X(s)N(s) + Y(s)D(s) = I, Re(s) > 0. (4) 

(N,D) is called a right coprime factorization (r.c.f.) for G. Left coprime factorizations 
(l.c.f.’s) are defined similarly. ( N,D ) is a l.c.f for G if G = D~ l N where N,D £ M( Hoo ) 
and there exists X,Y G M( “Hoo ) with 

N(s)X(s) + D(s)Y(s) = J, Re(s) > 0. (5) 

Every system which is described by a system of linear time-invariant ordinary differential 
equations has both a left- and a right-coprime factorization. (This is a consequence of the fact 
that the transfer functions of such plants are composed of rational functions.) Furthermore, 
the set of all stabilizing controllers for such plants may be described in terms of th eYoula 
parameterization ; a controller H externally stabilizes G if and only if it can be written 

H = (Y - RN)~ l (X + RD) y \Y -RN\f;0,R£ M(Hoo) (6) 

where X,Y,N,D are as defined in (4), (5). The set of all stabilizing controllers for a given 
system G are parameterized by R , as R ranges over all stable systems. In other words, 

S(G) = {{Y - RNy'iX + RD), \Y - RN\ ± 0,R £ M( Hoo )}• (7) 

The above formulation (7) is in terms of left coprime factors of the stabilizing con- 
trollers. The same family, S(G) may also be written in terms of right coprime factors: 

S(G) = {{X + DQ)(Y - NQ)~\ \Y-NQ\^0,Q£ M{ Hoo )}• (8) 

Unfortunately, more general systems, which do not have rational transfer functions, do not 
necessarily have either a left- or a right-coprime factorization. However, any system with a 
transfer function which can be written as a fraction ND~ l , TV, D £ M( Hoo ), is stabilizable 
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if and only if it has right (and left) coprime factorizations [16] . For systems of the type (1), 
the theory is more complete. First, some definitions are required. 

Definition 1.1: The (70-semigroup T(t) is stable if there exist constants M and 
a > 0 such that ||T(i)|| < Me~ at for all t > 0. 

Definition 1.2 A semigroup control system (A, B , C) is said to be internally stable if 
the semigroup generated by A is stable according to Definition 1.1 . 

Definition 1.3: The pair (A, 2?) is stabilizable if there exists a bounded linear oper- 
ator K : X — » R m such that A — BK generates a stable semigroup. 

Definition 1.4: The pair (A, G) is detectable if there exists a bounded linear operator 
F : R p — > X such that A — FC generates a stable semigroup. 

Definition 1.5: The system (A, B,C) is jointly stabilizable/ detectable if (A, B) is 
stabilizable and (A, C ) is detectable. 

The extensions to cr-stable, cr-stabilizable etc. are straightforward. 

Jointly stabilizable/detectable systems (A, 23,(7) have both left and right coprime 
factorizations [5]. Let K be such that A — BK generates a stable semigroup and F such 
that A — FC generates a stable semigroup. Define 

N(s) = CR(s, A — BK)B, (9) 

D(s) = / - KR(s, A - BK)B, (10) 

X(s) = KR(s, A - FC)F, 

and 

Y(s) = lAKR{s,A-FC)B , 

XN + YD - I. 

(N, D) is a r.c.f for the system i.e. CR(s ; A)B = AD -1 [10] and the Youla parameterization 
describes all stabilizing controllers for these systems. A l.c.f. can be defined similarly, using 
an operator F such that A — FC generates a stable semigroup . Define 

N(s) = CR(s, A - FC)B 

b{s) = I -CR{s,A- FC)F. 

Then, G(s ) = £>- 1 (s)iV(s) and 

NX + DY — I 
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for some X, Y G M( 7i ^ ). 

It is easy to show, using the Hille-Yosida Theorem, that internally stable control 
systems (1) are externally stable. We are also assured that external stability implies internal 
exponential stability: 

Theorem 2.1 [10] A jointly stabilizable/ detectable semigroup control system is internally 
stable if and only if it is externally stable. 

This equivalence justifies the use of controller design techniques based on system 
input/output behaviour for infinite- dimensional systems of the form (1). The discussion 
above is summarized in the following theorem. 

Theorem 2.2 Assume ( A,B,C ) is jointly stabilizable/ detectable and let G be the transfer 
function of the system. We say that a controller with transfer function H stabilizes G if the 
closed loop matrix A (G,H) G M(H 0 o ). Then H stabilizes the system, if and only if , its 
transfer function can be written using the Youla parameterization (7) (or (8 ) ). 

Also, if the state-space representations of the systems G and if are each stabilizable 
and detectable, then H stabilizes G , according to the definition at the beginning of this 
section, if and only if the closed loop is also internally stable. 

Theorem 2.3 [10] Assume ( A,B,C ) is jointly stabilizable/detectable and that H is a con- 
troller with a jointly stabilizable/detectable realization ( A c ,B c ,C c ). The closed loop system 
is input/output stable if and only if it is internally stable. 

The following is by now well-known, and several proofs exist. One version, which uses 
coprime factorizations, and which first appeared in [13] is given here. 

Theorem 2.4 Every jointly stabilizable/detectable system (1) is stabilizable by a finite- 
dimensional controller. 

Proof: Both the stable factors N, D have well-defined limits at in fini ty- 
lim a(N(s)) ~ 0, lim a(D(s )) = 1. 

p|- oo V V 7/ ’ p|— oo v v ” 

ile(s)>0 Re(s)>0 

It follows (Mergelyan’s Theorem) that each factor can be approximated by a rational element 
of M{ He*, ): for any e > 0 we can find rational N n ,D n G M( Tt^ ), with 

II m - WIL < II D{s) - D n (s) IL < e. 
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For sufficiently small e, ( N n ,D n ) is also a coprime pair. Let (X, Y) be a l.c.f of any 
finite-dimensional controller which stabilizes N n D n ~ l with XN n + YD n = I. Then, if e 
is small enough, Y~ l X also stabilizes the infinite-dimensional system since XN + YD = U 
where U has an inverse in M( Hoo )• Rewriting, we have X(NU~ l ) + Y(DU~ l ) = I. Since 
(iW7 -1 , DU' 1 ) is a r.c.f. for the infinite-dimensional system (1) , it follows that every jointly 
stabilizable/detectable control system (1) is stabilizable by a finite-dimensional controller. 
□ 


3 Approximation Scheme 

In this paper we are interested in the case where a finite-dimensional controller is designed 
using a numerical approximation to the original system. 

Suppose we have a sequence of finite-dimensional subspaces X n C X where the norm 
on X n is that inherited from X. Define P n x as the orthogonal projection of x € X onto the 
finite- dimensional subspace X n . We assume the following: 

(Al) The projection operators P n converge strongly to the identity on the Hilbert space X. 
That is, for all x € X 

lim II P n x — a? || =0. 

n—KX> n 

For each X n , the approximating system is (A n , B n ,C n )- 

B n := P n B, G n := C\x n , 

and A n is an approximation to A which satisfies the two assumptions (A2) and (A3) below. 
Note that the operators A n ,B n ,C n and the semigroup T n (t) generated by A n are operators 
on X n . 

A core C of a closed operator A is a linear space contained in the domain of A with 
the property that the set of elements (x,Ax), x <E C is dense in the graph Q(A) of the 
operator A ([11], pg- 166). 

(A2) We assume that there exists a core C for A such that 



for all x € C. 


Such an approximation scheme is said to be consistent. 


( 11 ) 
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(A3) We will further assume that the semigroups T n (t) generated by A n are uniformly 

bounded, that is, there exist real numbers N, M > 1 and k such that 

||r„(t)|| < Me kt for all n> N. (12) 

(Uniform boundedness of the approximate semigroups is generally referred to as “sta- 
bility” in the numerical analysis literature.) 

Consistency and uniform boundedness are sufficient for convergence of the approxi- 
mation [11, 14] i.e., for all e > 0,t > 0 and for all x G X there exists N such that 

II PnT{r)x - T n {r)P n x\\ < e for all r G [0 ,t] and n > N. ( 13 ) 

Assumptions (Al)-(A3) are satisfied by typical approximation methods. They are 
sufficient to ensure that the open loop response of the systems G n approximate the response of 
G. (Al) is of course redundant in that it is implied by (13). Assumption (A2) may be replaced 
by a requirement of strong convergence of the resolvants. (Trotter-Kato Approximation 
Theorem eg. [14]). However, these assumptions are not sufficient to ensure convergence 
of the closed loop response. An example of a scheme which approximates the open loop 
response of a system but is not satisfactory as a basis for controller design is given in [4]. 

The appropriate topology in which to establish convergence of the approximations is 
the graph topology. The importance of the graph topology in controller design is due to the 
following result: A family of plants G n can be robustly stabilized by a compensator H which 
stabilizes a nominal plant G if and only if G n converges to G in the graph topology. Fur- 
thermore, in this case, the closed loop response of the feedback pair A (G n ,H) converges to 
A(G, II). In order to obtain convergence of the closed loop response A(G'„, H n ) to A(G', H) 
for some H G S'(G'), Il n must also converge to H in the graph topology. The results given 
in this paper can be shown directly using the coprime factorization theory outlined in the 
previous section. The graph topology is not directly discussed further in this paper, although 
it is the topology underlying the convergence results. The interested reader is referred to [19] 
for an introduction, [20] for extensions to general algebras, [12] for its application to approx- 
imation of semigroup control systems, and [6] for robust stabilizablity under perturbations 
in the coprime factors. 

Suppose that the approximation scheme also satisfies assumptions of uniform stabil- 
izability, and uniform detectability: 

(A4) If the original system is stabilizable, then the approximations are uniformly stabi- 
lizable: there exists a uniformly bounded sequence of operators {K n } such that for 
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sufficiently large iV, the semigroups S n (t) generated by A n - B n K n are uniformly 
bounded by Me~ at for some M > 0, a > 0 and all n > N. 

(A5) If the original system is detectable, then the approximations are uniformly detectable: 
there exists a uniformly bounded sequence of operators {F rl } such that for sufficiently 
large N, the semigroups S n {t) generated by. A n - F n C n are uniformly bounded by 
Me~ at for some M > 0, a > 0 and all n > N. 

Common approximation schemes for several important classes of systems which sat- 
isfy these assumptions are given below. Not only do both approximation schemes discussed 
below satisfy assumptions (A1)-(A5), there exists a sequence of operators I< n which sat- 
isfy assumption (A4) and which converge strongly to an operator I< 6 B(X,U). Similar 
convergence exists for a sequence of operators F n which satisfy assumption (A5). 


3,1 Hereditary Systems 

Consider the delay functional differential equation 

x(t) = j drj(r)x(t + r) + Bu(t), t> 0 (14) 


y(t) = Cx(t), 


x(0) = x 0 , x(r) = (j>(T), -h < t < 0 

where x{t) R n ,y(t) € R p ,u(t) <E R m and B and C are matrices of appropriate dimension. 
Also, tj(t) is function of bounded variation taking values in R nXn with r/(r) = 0 and tj(t) 
is left continuous for -h < t < 0. Defining the state-space X = R n x C 2 (~h,0] R n ), the 
hereditary system (14) can be formulated as a control system of the form (1) [2]. 

Define the finite- dimensional subspaces of X to be 


X» = k«eX; Hr )=*>, -J^h < T < Fi-Jlh, j = 1 ATwherezj6.fr 


The finite-dimensional Galerkin approximation to (14) derived using these subspaces satisfies 
assumptions (A1)-(A3) [2, 15]. This scheme is known as the averaging approximation to (14). 

In [8] it is shown that this scheme is uniformly stabilizable and detectable: assump- 
tions (A4) and (A5) are satisfied . It is also shown that the stronger version of (A4) and 
(A5) holds. For a bounded operator I< such that A - BI< generates a stable semigroup, we 
can choose K n = KP n , K n satisfies (A4) and lim n -Mx> || K n - K\\ = 0. Similarly, if (14) is 
detectable, we can choose a sequence F n -> F where A- FC generates a stable semigroup 
and F n is a sequence which satisfies (A5). 
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3.2 Sectorial Operators 

Define A through a continuous sesquilinear form a :V xV h^C 


(~M, V>) = </>), <j> € V(A ), tp ev. (15) 

where (•, •) indicates the inner product on X, V is a Hilbert space , V(A) C V and V is 
densely and continuously in X. In order to avoid confusion with the norm on X , the norm 
on V will be indicated by || || F . Identify X with its dual so that V^X = X'^V'. We 

assume that in addition to (15), a(, ) satisfies Harding’s inequality: there exists k,c> 0 such 
that 

a(u,u) + k(u,u) > c|| u\\ v \ ( 16 ) 

The inequalities (15) and (16) guarantee that A generates an analytic semigroup with bound 


II 2X011 <e fc< - 

Further details may be found in [17]. Assume that the approximating subspaces X n satisfy 
a V- ap p r oxi in a. t i on property: for all x 6 V there exists a sequence x n 6 X n with 

Jiis, II - x llv = o. 

It is shown in [12] that this class of semigroup control systems, with a sequence of ap- 
proximating subspaces which satisfy assumption (Hi), leads to a sequence of approximating 
control systems which satisfy assumptions (Al)-(A5) of the previous section. Furthermore, 
as in section (3.1) we can choose K n ,F n so that K n I{ and F n -» F. 

Note that it is only required that projections onto X n converge in the V-norm. A 
similar result in [8] requires an inverse approximation property: 

II R ( s > A ) z ~ x Wv < e i( n )|| 4x 
II A * ~ x Wv — e 2 ( n )ll z \\x 

where ei(iV), e 2 (AT) -+ 0 as JV — > oo. These conditions are stronger, and more difficult 
to verify than (Hi). Furthermore, the above conditions can only be satisfied if R(s\ A) is 

compact. Problems such as control of structural vibrations with Kelvin- Voigt damping are 
thus excluded. 
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4 Controller Convergence 


It will now be shown that if both (A4) and (A5) hold in addition to the usual (Al)-(A3), 
then for large enough model order, the finite- dimensional controllers stabilize the semigroup 
control system (1). If a stronger version of (A4) and (A5) holds, then the controllers con- 
verge uniformly ( in M( H 0 0 ) to an infinite-dimensional controller, and the closed loop 
performance obtained with the finite-dimensional controllers converges to the closed loop 
performance obtained with the infinite-dimensional controller. 

Theorem 4.1 Assume the semigroup control system (1) is jointly stabilizable/ detectable 
and that the approximations used in design of the finite-dimensional controllers (3) satisfy 
assumptions (A1)-(A5). Denote the transfer functions of the controllers by H n and the 
transfer function of (1) by G. Assume, in addition, that the feedback operators K n converge 
strongly to an operator K 6 B(X,U): 

lim K n P n x = Kx , x € X. 

n— ► oo 

Then, for sufficiently large n, H n stabilizes the infinite-dimensional system G. For such n, 
the operators 

generate stable semigroups. 

Proof: Define 

X n (s) = K n R(s, A n - F n C n )F n 

and 

Y^(s) = I + K n R(s , A n — F n C n )B n . 

It is well-known (eg. [19]) that (X n ,Y n ) is a l.c.f. for the finite-dimensional controller H n , 
and that 

X n N n -f Y n D n = / 

where ( N n ,D n ) is the r.c.f of G n defined by (9) and (10) using the stabilizing feedback 
operator K n . 

It is clear from the assumptions that 

ii [*„ niL < m 

for some constant M. Let ( N,D ) be the r.c.f. of G defined by (9) and (10) using the stabil- 
izing feedback controller K. It was shown in [12] that lim^oo N n = N and lim^oo D n = D 


A ~BK n 

F n C A n — B n K n — F n C n 


(17) 
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and so, for sufficiently large n, say n > n 0 , 


|| [X n Yn] 


N-N n 

D-D n 


loo < 


n > n 0 . 


It follows that 


U n = X n N + Y n D 


has an inverse in M( Hoo ) for n > n Q . 

The closed loop matrix A (G,H n ) can be written in terms of G = ND~ X and H n 
Y~ x X n as 

I - NU~ l X n -NU-'Yn 


A (G,H n ) = 


DU~ l X n DU~ l Y n 


Since A (G,H n ) £ M(H oo ) the closed loop system is, by definition input/output stable. It 
is trivial to show that A on is the infinitesimal generator for the closed loop system. Internal 
stability of the closed loop follows from the fact that the controller is stable and Theorem 
2.3. Hence, A on generates an exponentially stable semigroup. Alternatively, it can be shown 
directly that the closed loop system is stabilizable and detectable and Theorem 2.1 applied. 
□ 


The dual situation, where K n is bounded and F n converges strongly to an operator 
F , can be shown identically, using convergence of a sequence of l.c.f’s of the systems G n 
and boundedness of r.c.f.’s of the compensators H n . We also have the following: 


Corollary 4.2 Let a be as in assumptions (A4)-(A5) and let a be any real number , 
a > a > 0. For n sufficiently large , the operator A on (17) generates an a-stable semigroup. 


Proof: As a consequence of assumptions (A4) and (A5) we can replace input/output stability 
and internal stability by input/output cr-stability and internal cr-stability in the previous 
theorem. The algebra Hoo is replaced by the shifted algebra Hoocr- n 

Theorem 4.3 Assume that the approximation scheme satisfies assumptions (A1)-(A5) and 
that the compensator design sequence is convergent: 


lim K n P n x = Kx , 

71-+O0 

for all x £ X and 

lim F n y = Fy 

n—>oo 

for all y £ R p . Let H denote the transfer function of the infinite- dimensional controller 
whose state-space realization is 

z(t) = (A - BK)z(t) + F(u(t) - Cz{t )) (18) 
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y(t) = Kz(t). 

The controller transfer functions H n converge in M( 7ioo ) ( he. uniformly in s) to H , 
Proof: We first show that l.c.f. ’s of the controller converge: 

lim || X - XnlL = lim II Y - KIL = 0. 

Define X n ,Y n as in Theorem 4.1 . Denoting the semigroup generated by A n — F n C n by 
S Fn (t )> the time-domain representations of these transfer functions are 

X n (t) := K n S Fn (t)F n , 

V n (t ) := S(t)I + K n S Fn {t)B n 
where 6(t) denotes the Dirac delta distribution. Define 

X(t) := I<S F (t)F, 

Y(t ) := S(t)I + I<S F (t)B. 

The proof is identical to that used to show convergence of coprime factors of the approxi- 
mating systems in [12]. Due to convergence of K n to K, and F n to F, we can show that, 

noo 

n^+oo j 0 II X n {t) ~~ X(t)\\ Rm +pdt — 0, 

roo 

„limjf o || Y.(t) - Y{t)\\ R ^A = 0 

where R m+p indicates some norm on operators from R m to R p . Convergence in this norm 
implies convergence of the Laplace transforms of these functions in M( Ti^ ). Each controller 
H n and the infinite-dimensional controller H are stable (elements of M( Ftoo ))• This implies 
that, Y n , Y are invertible in M{Hoo) and since lim n _>oo || Y n — YW^ = 0, 

lim ||yn" 1 -y- 1 || 00 = o 

and H n converges to H in norm. (This reflects the fact that the uniform topology is the 
restriction of the graph topology to stable systems.) Since (X, y) is a l.c.f. for the infinite- 
dimensional controller H , the result follows. □ 

Note that the controllers actually converge in the stronger norm: 

POO 

Km J || H n {t) - H{t)\\ Rm+p e^dt = 0, 
where a > 7 > 0 (a as in (A4) and (A5)) although this is not used here. 
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Corollary 4.4 The closed loop operator A(G n , H n ) converges in M{ Hoo ) to A (G,H). 
We can write the closed loop operator A(G, H ) in terms of stable coprime factors 

A(G,B) 


as 


and similarly, 


A (G n ,H n ) 


' I -NX 

-NY ' 

DX 

BY 

I - N n X n 

—N n Y n 

D n X n 

D n Y n 


The result now follows from convergence of the coprime factors. □ 

Thus we have shown uniform convergence of the controllers to an infinite-dimensional 
controller, and similar convergence of the closed-loop systems. 


4.1 LQG Optimal Control 

Commonly, a compensator of the form (18) is designed in two steps, through construc- 
tion of an optimal regulator and an optimal detector. Assume that (1) is jointly stabiliz- 
able/ detectable, so that construction of an internally stable closed loop is possible. 

If Q is a self-adjoint positive semi-definite operator and R is self-adjoint positive 
definite, then the optimal solution of 

POO 

J(u) = / (Qx{t), x{t)) x + (Ru(t), u(t))dt 

J 0 

is given by the feedback u(t ) = -Kx(t) where if := B* R^TL and II is the unique, non- 
negative, self-adjoint solution of the algebraic Riccati equation, 

( A* II + IIA — UBR -1 B* II + Q)z = 0, for all z € D(A). (19) 

The operator A — BK generates a stable semigroup on X. 

Similarly, we may obtain the observer gain F by solving the dual problem 

AS + S A* - E G* R-'CT, + Q 0 = 0 (20) 

where Q 0 is self-adjoint positive semi-definite and R 0 is self-adjoint positive definite. Setting 
F = % C* R~ x where S is the unique, nonnegative, self-adjoint solution to the above Riccati 
equation, A — FC generates a stable semigroup on X. 

Since the above infinite-dimensional Riccati equations cannot be solved exactly, ap- 
proximate solutions are calculated by using the approximations A n ,B n ,C n . Thus, if the ap- 
proximations are stabilizable, we obtain I< n = R~ l B* H n where II n solves, for Q n = P n QP n , 

A n II n + Tl n A n II n B n R * B* II n ■+■ Q n = 0. (21) 
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Similarly, if the approximations are detectable, the operators F n — E n C* R a 1 are found by 
solving, for Q on = P n Q 0 Pn > 

A n Yj n 4- A* — E n C 1 * R 0 1 (7 n S n + — 0* (^) 

The theory on convergence of solutions of approximating Riccati equations is fairly 
complete, and we can state the following which is a consequence of [3] and the theorems in the 
previous section. Convergence of the solutions requires that the adjoint semigroups S n (t) 
converge strongly, uniformly on bounded intervals, to S*(t). This assumption, together with 
(A1)-(A5) implies the stronger versions of (A4) and (A5) required for Theorems 4. 1-4.3. 

Theorem 4.5 Assume that the approximation scheme satisfies the usual (Al)-(AS) and that 
it is uniformly stabilizable (Af) and uniformly detectable (A5). Assume also, that for every 
x 6 X, the adjoint semigroups S* ( t ) converge strongly, uniformly on bounded intervals, to 

S* (t): 

lim sup || S* ( t)P n x — P n S* (t)a:|| =0, T > 0. (23) 

n-+i:>0 0<t<T 

Then: 

1. The controller (3) obtained by solving the finite-dimensional ARE’s (21) and (22) to 
obtain K n ,F n stabilizes the infinite- dimensional system for all sufficiently large n. 

2. This closed loop system is internally a -stable, for any a < a. 

3. The controller transfer functions converge uniformly (i.e. in M( Hoo )) to the trans- 
fer function of the optimal regulator/observer (18) obtained by solving the infinite- 
dimensional ARE’s. 

f. The closed loop operator A(G n , H n ) converges uniformly to A (G,H). 

Proof: It has been shown [3] that these assumptions are sufficient to guarantee that 
lim n _oo K n x = Kx and that lim^ FnV = Fy and that these feedback operators yield 
uniformly stable semigroups : (A4) and (A5) are satisfied. The conclusions now follow 
immediately from Theorems 4. 1-4.3 . □ 
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5 Conclusions 


A similar result to Theorem 3.1 and Corollary 3.2 has been shown in [8]. However, the use of 
coprime factorizations leads to a considerably shorter proof. Also, the results in [8] depend 
on several assumptions that are not required here. Unlike [8], the spectrum determined 
growth assumption not required at any point. Furthermore, both classes of approximation 
schemes discussed in [8] have the property that the resolvants of the approximations converge 
in norm to the original resol vant: 

Jim || 4)||=0 

and this property is used in the proofs. This restricts the class of systems under consider- 
ation to those with compact resolvants, and the approximation scheme to those with this 
convergence property. As demonstrated here, these assumptions are not necessary. 

The results here showing uniform convergence of a sequence of controllers, and of the 
corresponding closed-loop systems, have not been previously obtained. Gibson and Adamian 
[7] study the problem of convergence of LQG controller design for flexible structures, when 
the controller is designed using approximations to the solution of the partial differential 
equation describing the vibrations. Numerical results in [7], for finite-element approximations 
to structural vibrations indicated convergence of the controller transfer functions, but the 
authors were unable to prove Theorem (4.3) in the general case. 

One important consequence of convergence of the controller transfer functions, is that 
it justifies the use of such order reduction methods as Hankel norm reductions, balanced trun- 
cations, and coprime factor reduction. The order of the approximation (2) could be increased 
until the controllers (3) have converged sufficiently in M( Woo )• The numerical results in 
[7] for control of structural vibrations indicate that at least in some cases, convergence is 
obtained fairly quickly. The controller order can then be reduced using a standard technique 
such as balanced realizations. 
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